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INTRODUCTION 


Boundary layer skin-friction drag on an aircraft accounts for 
approximately A0% of the total vehicle drag. Therefore, any substantial 
reduction in skin-friction drag would lead to great saving in fuel. Skin- 
friction drag is induced at both laminar and turbulent conditions. Laminar 
flow (low skin-friction) is limited to surfaces near the leading edge, 
v^hile turbulent flow (high skin-friction) is dominant on the rest of the 
vehicle. Thus, any effective reduction in skin-friction drag should come 
from that part due to turbulent flow, and it has been suggested that this 
may be accomplished by use of compliant surfaces. 

Skin-friction drag reduction using compliant surface is a unique 
concept. Its main advantage is simplicity, no slots, ducts, or internal 
equipment of any nature are used. It requires only a surface with certain 
characteristics which would interact favorably with the flow, Kramer [1] 
is credited with the original idea of drag reduction by using compliant 
surfaces, based on his observations of dolphins swimming in water. 

Numerous experimental and theoretical studies have been undertaken to 
explain the phenomenon of drag reduction by compliant surfaces. A list of 
references and discussion of such results are given by Ash, et al . [2,3,4]. 
So far, there is no theory which can explain satisfactorily such a pheno- 
menon. However, some mechanisms have been postulated as possible explana- 
tions. One postulate states that the motion of the compliant surface 
delays the transition from laminar to turbulent boundary layer. Another 
assumes that the motion of the compliant surface alters the structure of 
the fully turbulent layer. The latter postulate appears to be more likely. 
Accordingly, Ash, et al. [3] have proposed two models for the alteration 


SUMMARY 


Some tentative compliant wall structures designed for possible skin 
friction drag reduction are investigated. Among the structural models 
considered is a ribbed membrane backed by polyurethane or PVS plastisol. 

This model is simplified as a beam placed on a viscoelastic foundation as 
well as on a set of evenly spaced supports. The total length of the beam 
may be either finite or infinite, and the supports may be either rigid or 
elastic. Another structural model considered is a membrane mounted over 
a series of pre-tensioned wires, also evenly spaced, and the entire mem- 
brane is backed by an air cavity. All these structural models belong to the 
general class of periodic structures for which a simple mathematical 
analysis is possible. The forcing pressure field is idealized as a frozen 
random pattern convected downstream at a characteristic velocity. The 
results are given in terms of the frequency response functions of the 
system, the spectral density of the structural motion, and the spectral 
density of the boundary- layer pressure including the effect of structural 
motion. These results are used in a parametric study of structural con- 
figurations capable of generating favorable wave lengths, wave amplitudes 
and wave speeds in the structural motion for potential drag reduction. 
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of the structure of the boundary layer. In the first model, the change of 
boundary-layer structure is attributed to the direct coupling interaction 
at the fluid-wall interface, and in the second model it is assumed that 
the surface radiates acoustical disturbances which alter the distribution 
of turbulent energy from the large-scale eddies to the smaller ones. 

In both models proposed in [3], desirable interactions of the surface 
with a turbulence field can be induced in one of three. possible ways; using 
an active wall, a passive resonant wall, or a flov/- triggered wall. 

The first type is one that is mechanically or electrically driven at 
prescribed frequency, amplitude, wave speed and shape. The objective is 
to obtain a forced narrow-band surface motion at a predetermined frequency 
and amplitude- 

The second type, a passive resonant compliant v/all , has been the most 
coimtonly-used experimental configuration in the past. In this case, the 
compliant wall is excited by the pressure field to establish a resonant 
mode. A favorable interaction would occur if this resonant mode would 
vibrate in a frequency that falls within the energy-containing spectral 
range of the turbulence. 

The flow- triggered wall (.or in the true sence of the word, an ideal 
compliant wall) is one that responds to all the turbulence "bursts", not 
only to a few "bursts" as in the other two types. Although this kind of 
a compliant surface may be the most effective in drag reduction, it is 
extremely difficult or even impossible to design such an ideal surface. 

Several possible structural configurations have been suggested for 
passive compliant surfaces: membrane surfaces, rigidly backed slabs, 

laminated structures, and periodic structures (structures which are 
composed of identical units repeated at equal distances). Of all these 
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configurations, periodic structures have received the least attention. 

In this investigation, therefore, the design of periodic structures as 
tentative compliant surfaces will be emphasized. In particular, two 
models of periodic structures will be investigated: a ribbed membrane 

backed by a soft material such as polyurethane or PVS plastisol, and a 
membrane backed by an air cavity. In the first model, the structure is 
assumed to be resting on a viscoelastic foundation which v;ould approximate 
the behavior of polyurethane or PVS plastisol mentioned above. Also, in 
this analysis, the ribbed two-dimensional membrane is replaced by a one- 
dimensional periodic beam on equi-spaced supports. The supports, which 
are simplified versions of ribs, provide rotational and transversal 
elastic restraints to the beam. 

In the second model, a number of circular solid wires are attached 
(glued) to a two-dimensional membrane. In the analysis, the wires are 
treated as pre- tensioned strings. 

In either model the structure can be of finite or infinite length. 

For the first model, both cases will be considered while only an 
infinitely long structure v/ill be considered for the second one. 

In order to achieve drag reduction, the proposed design should 
possess certain dynamic characteristics which will interact favorably 
with the flow. In other words, the motion of the compliant surface when 
excited by a turbulence field should have the wave shape, propagation 
speed, wave length, and amplitude which would make such interaction 
possible. In the absence of a theory which can predict such characteris- 
tics, wind tunnel measurements have been the only available source of 
information. Historically, experimental investigations have been carried 
out at low subsonic speeds. This is probably because of the fact that it 



is easier to control the experiment at such speeds. Accordingly, Ash, 
et al. [3,4] have done some experimental investigations at a speed range 
of 50-150 ft/sec to verify some of the previous experiments on one hand, 
and to try to establish a compliant drag reduction theory on the other hand. 
Based on their wind tunnel measurements, they have suggested that the 
desired motion of a possibly-successful compliant surface should have 
the following characteristics- (Private communication from R. ,L. Ash) 

-3 

(1) the maximum amplitude of the vibration should be (1-3) x 10 
inches 

(2) the resonant frequency should not exceed 1500 Hz 

(3) the wave length of the surface motion should not exceed 
.1 inches. 

In this analysis, the above characteristics will act as criteria 
for the design of a tentative compliant surface. However, one should 
keep in mind that the actual goal of compliant drag reduction is to 
reduce drag not only at such low speeds, but most importantly at high 
subsonic and supersonic speeds which represent a range within which 
most vehicles fly. 

The response of finite beams to a convected pressure field has 
been analyzed in the past by calculating the principle modes’ of the 
system, and by studying the forced motion in each significant mode [5-8]. 
Alternatively, transfer matrix techniques have been developed by Lin, et 
al. [9-12] to calculate the response directly without determination of 
normal modes. In the case of periodic structures, accurate computation 
of normal modes is extremely difficult. Mead and associates [13-17] 
have used another approach, namely the wave propagation approach which 
is originally due to Brillouin [18]. This approach works well and with 
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great simplicity for periodic structures of infinite length. It can also 
be applied to finite periodic structures but the procedure becomes more 
involved. , 

In this investigation, the wave propagation approach also will be 
used to obtain a closed-form solution for the structural response. 

However, the analysis differs from that of Mead in two respects: first, 

it takes into consideration the effect of the foundation (or the air gap) 
on the structural response, and secondly ^it incorporates the effect of the 
interaction between the fluid and the motion of the structure. Finally, 
the result of a parametric study aiming at an optimal design which has 
the desired surface motion will be presented. 
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I. THE EXCITING PRESSURE FIELD 


The structure is assumed to be excited by a turbulent pressure field 
stemming from a fully turbulent boundary layer. The pressure field is 
random in nature and can be described only by statistical quantities. 

From the standpoint of the structural response, the simplest mathematical 
model for such a field is one that is statistically homogeneous in space 
and convected in a given direction as a frozen pattern. The frozen pattern 
assumption is known as Taylor's hypothesis. Since experimental measurements 
of the cross-correlation or cross-spectrum of a boundary layer turbulence 
generally show some degree of decay in addition to convection, Lin and 
Maekawa [19] recently have proposed a superposition scheme where a decaying 
turbulence can be constructed from infinitely many frozen-pattern components. 
Although the new model might represent a real turbulence more accurately, 
Taylor's hypothesis will be used in this investigation for computational 
savings. 

Therefore, the turbulent pressure field will be treated as a random 
function of x - U t and y where U_ is the convection speed, t is time, 
and the flow is assumed to be convected in the positive x-direction. 

Under this assumption, the pressure over a panel of width b can be 
expressed as a Fourier-Stieljes integral 

GO 

P(x - U^t.y) = I e^(“t - kx) 

-oo 

vihere 

(0 = frequency in radians per second 
k = wave number, {k - •^) . 

It can be shown that since P is statistically homogeneous 
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E{dF(k-|,y^) dF*(l<2,y2)> = Sp(kpn) 5(k^ - k23 dk^dkg 


where n = y-j - y^, £{♦} represents the ensemble average, an asterisk 
denotes the complex conjugate and Sp(k,ri) is the wave-number cross-spectrum 
referred to a coordinate frame moving with velocity To simplify further 
our analysis, we assume that can be expressed in a separable form as 

Sp(k^,n) = 


where g(n) is the cross-correlation of turbulence in the y-direction. In 
some cases the functional form of g can be determined analytically. For 
instance, in the case of homogeneous and isotropic turbulence, g{n) is 
given by Lin [20] 


9 ( 0 ) = k 


1 


2 

n y fari\ 
2 *^2^ L ^ 


a 



3K, 


,an’v <in y. /aru 
V ~ L L ^ 


an 


raTii 


2 h 

where L is the scale of turbulence, a = [1 + (Lk-j) ] , and K^(*) is the 
Bessel function of the second kind (Although the above expression for g(n) 
is derived for the cross-spectrum of turbulence velocity, it is reasonable 
to use the same form for the pressure using a different value for L). 

However, based on experimental data for homogeneous boundary layer turbulence, 
Mastrello [21] expresses g(n) as 


g(n) = e-l’il/® 


( 1 . 2 ) 


where B is a characteristic length proportional to the boundary layer 
thickens 6. In what follows we will use Mastrello's expression, Eq. (1.2). 
Expanding Eq. (1.2) in a Fourier series. 
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(1.3) 


gb)= E 


n=-“ 


The Fourier coefficients f^ 's are computed from 


^ l'’ e-’"’'!'/” dr, 

-b 


1 

b 


cos n) dri 


(1/B) 


^ {1 - (-rf 


^ 0/3)^ + ® 


(1.4) 


Mow, the separability assumption of Sp is equivalent to replacing 


Eq. (1.1) by 


P(x - U^t.y) = G(y) 


f dF(k) 


(1.5) 


with 


E(G(y^) G*(y2)} = g(n) 


If G(y) is also expressed as a Fourier series 


oo <Xl 

G(y) = I f e 

n=-oo 


irniy/b 


( 1 . 6 ) 


then by use of the orthogonality of Fourier coefficients we can show that 

,2 
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In the case of a one-dimensional structure, the variation of a 
pressure field in the y-direction becomes immaterial. Therefore, the 
pressure P may be expressed as 

CO 

P(x - U^.t) = I dF(k) (1.7) 

»oa 

It can be seen from Eq, (1.5) or Eq. (1.7) that if a linear analysis 
is valid, the solution to the present random vibration problem can be 
constructed from the solution of a fundamental problem where the 
excitation is 

(a) a unit plane wave, for the one-dimensional case, 

(b) a y-direction modulated unit plane wave, ginTry/b^ 

for the two-dimensional case. 

It seems proper to call such a fundamental solution a frequency response 
function. 

In the following sections, we will be mostly concerned with obtaining 
the frequency response functions for different structural configurations 
corresponding to a unit excitation of the type either (a) or (b) mentioned 
above. In each case, the frequency response function can then be used to 
compute the spectral density of the structural response. The spectral 
analysis of the radiated pressure induced by the structural motion also 
will be presented. 
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II. AN INFINITE BEAM RESTING ON A VISCOELASTIC FOUNDATION 


The configuration of such a beam is shown in Fig. 1 (a and b). 

Fig. 1-a represents the case when the supports are considered to be 
infinitely rigid in the transverse direction, while Fig. 1-b represents 
the case when these supports are transversely elastic. As can be seen 
from Fig. 1, the structure is composed of an infinite number of bays, 
each of length A. Also, it is to be noted that the representation of the 
viscoelastic foundation as shown in Fig. 1 is just one of several possible 
representations [27-30]. The inclusion of the spring k-j is to accommodate 
some viscoelastic materials which deform instantaneously under loading. 

The values of the spring constants k-j , V .2 and the damping coefficient 
depend on the specific chosen material. 

The equation of motion of such a beam at any location not immediately 
over a support is given by 


where 


D ^ M + F = P + P, 


z=0 


( 2 . 1 ) 


P = the turbulence pressure field neglecting the effect of the beam 
motion 

P^ = additional pressure generated by the beam motion 
F = foundation reaction 
W = transverse displacement of the beam 
mjj = mass of the beam per unit length 
D = flexural rigidity of the beam. 


Since we are interested, at this stage, in obtaining the frequency response 
function only, a unit excitation will be used instead of the 







random pressure P. Also, if the periodic beam considered is actually a 
strip of periodically supported plate, the more appropriate value for D 
to be used in the computations would be the bending rigidity of the plate. 
By the same token, mj^ and F would be mass and force per unit area, 
respectively. 

The reaction of the foundation, F, is obtained as follows: 

From Fig. 2, one finds that 

W = + W 2 

where 

W-j = relative displacement of point a w.r.t. point b 
W 2 = relative displacement of point b w.r.t. point c 

(which is the absolute value of the displacement of point b) 

F = (2.2) 

F - ^2W2 + Hy W 2 (2.3) 

To obtain an expression relating the reaction F and the total deflection W, 

differentiate both sides of Eq. (2.2) 

F = k^ W , (2.4) 

multiply Eq. (2.2) by k 2 » Eq. (2.4) by n^, Eq. (2.3) by k^ and add, 

(k, + k2)F + ny F = kik2(W^ + W 2 ) + k^ riy (W, + Wg) 

or (k.j + kg)F + Hy F = k.jk 2 W + k^j W (2.5) 

.Eq. (2.5) can be written in a more suitable form (for future discussions 
concerning the effect of the spring k 2 ) by dividing both sides by k 2 

ki rj . k-j • 

(1 + ^)F + F = k,W + W (2.6) 

^2 ^2 * ^2 
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Fig. 2 Viscoelastic Foundation 
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Combining Eqs. (2.1) and (2.6), one obtains 


k, n„ 

{ 0 + Hj^) + ^ 


3 } (D ^ + m. ii) + k,W + ^ ft 

3x^ D I Kg V 


= ^ ) « 2 . 7 ) 

^ ^ z=0 

\ 

For an Infinitely long periodic beam, the wave propagation approach 
due to Brillouin and Mead can be applied ffiost conveniently. Rewriting the 
exciting unit plane wave in the form, 

jKut - kx) , ji(ut - Uq-X/I) (2 g, 


where Uq = kA 1s seen to be the phase difference of the exciting wave at 
two points separated by a distance £ along the beam. However, the structural 
response is not restricted to this phase difference although it must be 
spatially periodic with period £.. This periodicity condition Is satisfied 
by expressing the structural response as follows: 

W(x,t) = Z A e^^“^ ■ (2.9) 

n*-“ 


with = Pq + 2mr, This series must be made to satisfy the B.C, 's at the 
supports. Also, It is to be observed that corresponding to different 
different wave components travel along the beam at difference speeds; 
namely 


r - 

n " ~ Pq + 2nir 


( 2 . 10 ) 


Each wave component can travel in either the positive or negative direction 
depending on the sign of y^. For small values of IPpl. the absolute value 
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of can be very large. The physical implication is that a subsonically 
convected excitation can generate supersonic waves in the structure. On 
the other hand, subsonic waves also can be generated by supersonically 
convected excitations since become smaller for large values of n. 

Another interesting relationship involves the wave lengths of the 
response components. If the usual definition of wave length is modified 
so that a negative-going wave has a negative wave length, then the wave 
length, is related to by 

= constant 

n n 

Thus, slower traveling waves have shorter wave lengths. 

Now, to solve Eq. (2.7) for the displacement W, the reactive pressure 
P-j has to be expressed in terms of W. The pressure which is induced by 
the structural motion on the upper side of the beam is governed by the 
convected wave equation 

(A . V 1-)%, = a? 7^ P, (2.11) 

where V = the flow velocity 
a-j = the speed of sound 

I 

2 

V = Laplace operator 

This induced pressure is subjected to the boundary condition 



z=0 

and the radiation condition that it propagates in the domain z < 0. 

The solution to Eq. (2.11), taking into account the above conditions, is 
given by Morse and Ingard [22]: 
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(2.13) 


“ faPl^l (Cn - ■ 

1 I n- CO / 9 " 

Z=0 /(c„ - V)2 - 

where p-j ® fluid density surrounding the beam (air) 

Cjj * the propagation speed of the nth wave component of the response 
= nth component of the displacement 
. A„ ' 

Some general comments about the expression for as given in Eq. (2.13) 
are in order: 

(i) When = V, the contribution of such wave component to is zero, 

(ii) When IC^ - V| = a^, P-| -»■ <» (shock wave effect) 

(iii) The nth component of P-j will give a damping effect whenever 
jCn - Vl > a^. This would occur: 

(a) If u is very high which makes very large (especially 
when is negative) 

(b) When V > a-j and is negative (or is positive but 
|C^ - Vl is still greater than a-|) 

(iv) The nth component of P-j will act as an added mass if the quantity 
under the square root sign in Eq. (2.13) is negative; i.e., 

|c„ - VI < a, 

(v) When V < a^ and to is small (C^ — component of P-j 

2 

will be approximately proportional to V . 

Determination of the Coefficients A^'s: 

Substituting Eq. (2.9) into Eq. (2.7), one obtains 

D^+CeW = P + Pi| (2.14) 
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where 


Ce = 


(1 + iuk^Tiyk2 

1 + ^ + '“V=2 


(2.15) 


E w„ 

n=-oo 


Up to this point, nothing has been said about the boundary conditions 
at the supports. In the following sections, two cases will be considered: 

A. The supports are assumed to be infinitely rigid in the transverse 
direction. 

B. The supports are assumed to have some elasticity ia that direction. 
A. Transversely Rigid Supports 

If the supports of a periodic beam are rigid, the displacements at 
these supports are zero. Hence, at the support x = 0, 

W(0) = I A^ = 0 (2.16) 

n=-oa 

or 

Ao=- I A„ (2.17) 

n=-oo 

n?^0 


Thus, the displacement, W, and the pressure, P.j, can be expressed as 
(the common factor e^“^ will be dropped from now on for convenience) 

W = I A^(e"'‘V/^ - (2.18) 

n=-oo 

n?«0 


\ = I lup^a^ 

n=-oo 

nj«0 






•) K 


(2.19) 
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where 


\ ° - =1 >V (2-20) 

The coefficients A^,s may be obtained by using the virtual work 
principle (i.e., when all the forces on and in the beam are in equilibrium, 
they should do no virtual work through any virtual displacement). 
Specifically, we assume a virtual displacement 

. ( 2 . 21 ) 

As usual in complex algebra, the conjugate of the virtual displacement, 
namely 

5* = 6Ajjj (6"*^/^ - (2.22) 


is used in calculating the virtual work. Due to the spatial periodicity of 
the forcing function, the structural response and the virtual displacement, 
the result of an integration over an entire structure composed of B periodic 
elements is equal to B times that of integrating over one element. Since 
the total virtual work done is zero, the virtual work done within each 
periodic element is also zero and it is sufficient to apply the virtual work 
principle to just one periodic unit. 

Now, the contribution of the 1st term in Eq. (2.14) to the virtual 
work is 




e 


-ipQx/ij 


6A^ dx 


“ B [(li^J.)'’ + (ug/Jl)'' 



n=-oo 

nj«0 


3*S-*6A 


m 


(2.23) 
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Note that in obtaining Eq. (2.23), the following relationships have been 
used: 


(1) fV’v^*' • 

dx = < 

;o 

\ 

; m j« n 

(2.24) 

0. 



; m = n 


(2) fV’v/^ 

dx = 

t 

; m jf 0 

(2.25) 

0 



; m = 0 



Similarly, for the other terms in Eq. (2.14) 


6 W 


( 2 ) 


= ^e r ^ A ^( e "’’’ V /^ - • 6 A ^ 

■'n n=-» 


- e'^^0^''^) dx 


nj^O 


(2.26) 


6W. 


= - e-‘''V/^ (e'''V/^' - e^ V/^) dx 


■I 

'0 

0 


; m = 0 


I • 6A_ ; m 0 
m 


(2.27) 


6F„ 


rl ~ 

= - I 1(op,a, I A„. 
L n=-“ 


,-iyj,x/a g-iyo^^/A 


0 " I “n *"0 

JA^ (e'V/^ - dx 


loip-ja-j 




Kijp-ja-j » 


m 


I ^ " 

n=-oo 

n?«0 


(2.28) 
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The virtual work done by the rotational spring, k^, at the support x = 0 
is 

6Wp = kre^(0)(6Aj'(0) = k^W(0)(6Aj'(0) 


where the prime, ', denotes differentiation w.r.t.x. 


= 


n=-oo I 

<■ 

- X 


6A. 


m 


n?£0 


2nir 2tmr 


6A. 


m 


(2.29) 


The total virtual work, 6W, is 


6W = 5W,,x + 6^r,o^ + 6W„ + + 6W^ = 0 


‘(1) '^"(2) 


1 


(2.30) 


Using Eq. (2.23) and Eqs. (2.26)-(2.30) , we obtain a set of simultaneous 
linear equations for the coefficients A^'s as follows: 


ft j. r •*' 

A„ Um + -pp 

m m e ’’ 


top, a, j4\ - . / 4 . . '“<=1=1 s'*' 

/ n?«0 \ ° • 


^r D z 5, ^n D 
n — 

nj«0 


(2.31) 


The structural damping of the beam and the supports can be accounted for 
by use of complex D and ky. 


D = D'(l + inj,) 
k^ = kv(l + in^) 
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where and are the loss factors for the beam and the support materials 
respectively. Also, introducing the nondimensional quantities 


r.2 2 sr 

“ = "?b“ : 


Pi = P^wa^ ^ 


Y = k; A/D' 


a = A^/D' 

3 = kg A^/D' 
Ti = ojTiy A^/D-' 

^ = wHy/kg 




(2.32) 


Eq. (2.31) can be simplified as 


<|){m) + (j)(0) I , A^ + ' I A^C(n. m)=S 

n=-oo n=-“> 


n/0 


n/0 


(2,33) 


where 


^(m) = 


+ 1 ^ + iij^y + a + in ^ 

4 . ^ 

^m ^ 


(1 + + i^)(l + injj) 


+ ini,) 


(2.34) 


1 + in. 


5(n,m) = Y • I + ~ i7| ' ^ * dmr) • (2mrr) 


(2.35) 


S = - 


A^/D' 

1 + ini. 


(2.36) 


In practical applications the number of terms in each infinite sum in (2.33) 
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must be truncated. The truncated version of (2.33) may bte written in a 
matrix form as 


where 


G(A} = {S} 


+ 4>(o) 

= Kronecher delta = 

t.j = -N. -(N-1). - H-1. N 


(2.37) 


(2.38) 


The choice of the number of simultaneous equations (»2N) depends on the 
desired accuracy in computing the Values of A^^'s. 

In the remainder of this section, two special cases will be discussed. 
Iitj the first case, the supports will be assumed to offer no rotational 
v'estraint to the beam, and in the second case the foundation will be 
assumed purely elastic. 

(1) No rotational restraint at the supports (i.e. , kp » 0): 

In this case Eq. (2.33) becomes 


A„ 4>(m) + <^(0) I A^ * S (2.39) 

n=-oo 

nj^O 

The solution of this equation is much simpler than that of Eq. (2.33) 
as to be explained below. 

From Eq. (2.39), one observes that 

Ajji (}»(m) = constant 

or A^ (|)(1) = A 2 <j>(2) * *(8) (2.40) 

Hence, combining Eqs. (2.39) and (2.40) 
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^ = (2.41) 

1 + (j)(0) Z l/(j)(n) <J)(0) Z {l/(j)(n)} 

n=-oo n=-oo 

njfO 

which readily gives 

\ (2-«) 

<J)(m)<|)(0) Z l/<^(n) 

n=-oo 


The number of terms required in Eq. (2.42) to obtain an accurate total 
deflection can be simply determined. First, we establish the dominant term 
in the series. Referring to this term as A^, compute the ratio 




4>(nH) 


4>(n) 


(2.43) 


Defining the truncation as -N £ n £ N, a reasonable criterion for choosing 
N is 

r^^ « 1 for jn| > N . 

(2) Elastic foundation ; 

As can be seen from Fig. 2, the foundation can be made elastic by choosing 
(a) kg = “ 
or (b) Tiy = 0 

or both. 

Case (a) 

If k2 then Eq. (2.15) becomes 

Cg = -a)^m|j + k-j (2.44) 

and Eq. (2.34) reduces to 
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I 







+ a 




m 




(2.45) 


while 5(n,m) and S remain unchanged. 

Case (b) 

When Hy = 0, 

Ce = irrVir 


and 


4>(ni) = - 


1 ■'2 ^ 


as 


. P] 

1 — L 


1 + 1 % 1 ■*■ 


G„ 1 + iriu 
m b 


(2.46) 


(2.47) 


Again ^(n,m) and S are unchanged. 

Combination of cases (1) and (2) (i.e., no rational restraint at the 
supports and with an elastic foundation) will, of course, lead to further 
simplifications. The coefficients A^'s are obtained most simply from 
Eq. (2.42) with c|)(m) given by either (2.45) or (2.47). 

Once the coefficients A^'s are determined, the displacement W and 
other related quantities can be readily computed. For instance, the first 
derivative (w.r.t. x), W' , is given by 

W'(x) = I A |-(i ^) e"''V''^ + i ^ e"''V^A (2.48) 

n=-oo J 

and 

CO 

w’(0) = I (2.49) 

n=-oo 

r 

Similarly, the second derivative is given by 

W"(x) = I A„((-i ^)^ e"''V/^-<(-i ^)^ eS^^/A (2.50) 

n=-oo V, 
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(2.51) 


W"(0) = f (- 

n=-oo 

Since we have substituted a unit plane wave for the excitation 
oressure P, the computed response as given in Eq. (2.18) (with the time 
variation factor e^“^ omitted) is actually the frequency response function 
H(x,k). Therefore, we can write 

H(x,k) = I A^(e"''V''^ - (2.52) 

n=-» V J 

where the A^'s are obtained by solving the matrix equation (2.37) or 
directly from Eq. (2.42). 

B. Transversely Elastic Supports 

Eventually if a periodic structure is used as a compliant surface, 
the supports are most likely ribs or stringers which have inertia and 
stiffness. In this section, we will assume that these supports can be 
replaced by springs characterized by a spring constant k^ in the transverse 
direction as shown in Fig. 1-b. The inertia of the supports will be 
neglected. An analysis similar to that of the previous section will apply 
here. Avoiding repetition of the previous section, the differences may be 
pointed out: 

(1) The displacement at the internal supports are no longer equal 
to zero and the displacement is given by 

W(x,t) = I A e^“^ (2.53) 

n=-oo 

(2) There will be an additional contribution to the virtual work. 

This comes from the elasticity of the supports. Following the 
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same procedure as in the previous section, the virtual work done 


by at X = 0 is 

alTt = kt W(0) = 6A„ • f A„ 


(2.54) 


n=-oo 


Then the relationship governing the coefficients A^^'s becumes 
A„ <|)(m) + r X + yy„ i = / -S ; m = 0 

«" n=-co " n=-co " I 0 ; m 0 


(2.55) 


3/nt 


where r = k^ /D* , y and S are the same as before, and the expression for 
(j)(m) is given by Eq. (2.34). 

The first and second derivatives of the displacement are given by 


W(x,t) 

W'(0,t) 

W"(x,t) 

W"(0,t) 


w 

1 


(2.56) 

n=-oo 

00 

A_(-iu_ /I) e’^^ 


1 

(2.57) 

n=-oo 

00 

1 

A /t)^ e’“* 

(2.58) 

n=-oo 

00 

1 

- (m„/a) 2 A„ e'-“* . 

(2.59) 


n=-oo 


and the frequency response function H(x,k) is given by 
H(x,k) = I A^ e''*V/^ 


(2.60) 


n=-oo 


Equations (2.48)-(2.52) and (2.56)-(2.60) will be used in the following 
sections. 
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It is worth noting that analyses of certain special cases discussed 
in the previous section, including the case of zero rotational constraints 
at the supports and the case of the elastic foundation, are also applicable 
here with the obvious modifications which correspond to the change in the 
B.C.'s at the supports. 
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III. FINITE PERIODIC BEAM ON A VISCOELASTIC FOUNDATION 

The configuration of such a beam is shown in Fig. 3. Fig. 3-a shows 
a beam with transversely rigid supports, while Fig. 3-b shows the same beam 
with transversely elastic supports. 

The response of a finite periodic beam to a unit excitation, 
can be viewed as consisting of two parts: 

(1) The response of an infinitely long periodic beam with the same 
periodic element. Such a response has already been discussed in the 
preceding sections. 

(2) The response corresponding to the reflections from the end 
supports of the finite beam. 

The waves reflected from the left-hand end support are positive-going waves, 
while the ones reflected from the right-hand end support are, naturally, 
negative-going waves. The reflected waves are free waves and can be viewed 
as the solution of the equation of motion without the excitation wave (i.e., 
by putting P = 0 in Eq. (2.14)). In other words, the motion is a free 
vibration of a semi-infinite periodic beam. It is important to note that 
the total response, which consists of parts (1) and (2) explained above, 
should satisfy the boundary conditions at the extreme ends of the finite 
beam. 

The equation of motion for free vibration of the beam is given by 

D^+CgW-PJ =0 (3.1) 

z=0 

The reactive pressure P-j depends on the displacement W. For a semi- 
infinite periodic beam, there is no known closed-form solution for P-j. 
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z 


Fig. 3-a Finite Periodic Beam {Transversely Rigid Supports) 



z 

Fig. 3-b Finite Periodic Beam (Transversely 
Elastic Interior Supports) 
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Therefore, an iterative procedure will be used. Since the effect of is 

expected to be small, we will initially neglect such an effect (i.e., put 

P-|j = 0 in Eq. (3.1) ). Next, determine based on that solution of W, 
z=0 

and a new displacement will be computed by solving Eq. (3.1) with previously 
determined P^ treated as an external force. This will become clearer as 
we proceed. 

Rewrite Eq. (3.1) (after putting P^ = 0) as 

^.RM=0 (3.2) 

3x 

- ‘^e 

where R = - -g- . The solution of Eq. (3.2) can be written as 

4 XpX 

W(x) = jAj^e ;0<x<t (3.3) 


in which the X^'s are the four complex roots of ^ and they are 



(3.4) 


where Re, Im denote the real and imaginary parts of a complex quantity, resp. 

Now, consider a positive-going wave which has a propagation constant p. 
The response produced by such a wave in a periodic structure has the 
essential property of which: 
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vector at station (r+1) = x vector at station (r) 


(3.5) 


where r and (r+1) are two stations on the periodic beam separated by a 
distance equal to the length of the periodic element, i, as shown in 
Fig. 4. The components of such a vector can be taken as any four 
independent linear combinations of the displacement and its first three 
derivatives w.r.t.x. 

The coefficients A^^'s of Eq. (3.3) will be determined by using the 
wave property explained above, and the boundary conditions at the extreme 
ends of the beam. In what follows, two cases will be examined: transversely 

rigid supports, and transversely elastic supports. 

A. Transversely Rigid Supports 

For transversely rigid supports, the displacements at both, ends of 
each periodic element (x = 0, Z) are zero. Therefore, Eq. (3.3) gives 

4 

I A. « 0 (3.6) 

j=l ^ 

4 X .S, 

and ^ A. e '^ = 0 (3.7) 

j=l J 

Using the wave property expressed in Eq. (3.5) to relate the slopes at 
X = 0 and X = £, 

W'(£)=e^W‘(0) (3.8) 


Substituting Eq. (3.3) into Eq. (3.8), we obtain 


X -J, y 4 

w . 


I A. X e = e I A X 
j=l j=l 


(3.9) 
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station r 


station r+1 


Fig. 4 A Single Periodic Element 

(Transversely Rigid Supports) 


Furthermore, one can observe (from Fig. 4) that 


W'(il) = D W'{1) 


Since M(x) = D W"(x), Eq. (3.10) becomes 

k„ 

W"(il) + W'(t) 


or 


,2 . '^r 


.1^ A. e 


X.£ 

J ~ 


e^ W"(0) 

A 9 

' I ^ j 

j=l ^ ^ 


(3.10) 


(3.11) 


Equations (3.6), (3.7), (3.9) and (3.11) form a set of simultaneous equations 
for the coefficients A^'s. However, these equations are not sufficient to 
determine the coefficients A^'s because the propagation constant p is still 
unknown and the problem will reduce to an eigenvalue problem as to be 
demonstrated below. 

Combining Eqs. (3.6) and (3.7), one obtains 


where 

Substituting Eqs. 
where 


Ag = e'^ {-A3e3 - A^e^} (3.12) 

X.l 

e. = e ^ - e ' ; j = 2.3,4 (3.13) 

J 

(3.12), (3.13) into Eqs. (3.9) and (3.11) and simplifying, 

G {A} = e^ Q {A} (3.14) 

{A}=(?3V (3.15) 

^”4J 


G = 


^ 3 "^ 2®2 ®3 




®4 


‘^ 4"^2 


(3.16) 
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^3“^2®2 ®3 

^4"^2®2^ ®4 

-1 

-1 

93-92^2 e4 

94-92^2 


X.s. X,A 

f . = X. e ’ - X^ e ' 

9-j “ 

di = (x2+-^Xi) e ^ 

S-j ~ d-j + i^2 ~ ®2 ^3 

$2 ® di + (d2 ” d^ ) &2 6/^ i * 1 (2y3)4 

Rewriting Eq. (3.14) in a more suitable form as 

C {A} = e^{A} 

where C = Q~^ G 


(3.17) 


(3.18) 


(3.19) 

(3.20) 


Equation (3.19) is a standard eigenvalue problem in which e^ is equal to the 
eigenvalues of the matrix C. These eigenvalues, denoted by A-j and A 2 , are 
pairwise reciprocal [15] and they are 

^ 1,2 i {^11 ^22 - /(Cii-C22)^ + 4 ^21 } 

where C^. is the (i,j)th element of the matrix C. The propagation constant 
y is given by 

H = ± 1o9e = + l°9g ^2 

= + iyj (3.22) 
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The real part of the propagation constant, represents the decay of 
the amplitude of the propagating wave while Uj represents its phase shift. 

We can see from Eq. (3.22) that two waves propagating in opposite directions 
will be generated in the structure. Also, the eigenvectors corresponding 
to y and -y can be determined. Of course, only the ratio between the two 
components A, and A^ of each eigenvector is unique. This ratio can be used 
in conjunction with Eqs. (3.6) and (3.121 to determine the coefficients A^‘s 
in terms of one of them (say 1, Ag/A-j, A^/A^ and A^/A-j). If we denote 


Aj/A 


1 



; for +ve going wave 

1 for -ve going wave , 
j=l-*-4 


then Eq.- (3^3) becomes- 


W(x) = < 


4 X.x 
w+ I at e J 
j=l ^ 

. 4 X.x 

w I a. e ^ 

L “ ^ 


+ve going wave 


ve going wave 


(3.23) 


(3.24) 


where w_^ and w_ are the amplitudes of the positive and negative going waves, 
respectively, and will be determined by using the BrC's at the extreme ends 
of the beam. 

Now, to determine the induced pressure P-j, we have to solve Eq. (2.11) 
subjected to the B«C given by Eq. (2.12) and to the radiation condition 
that it propagates in the negative-z domain. For a positive going wave 
with y as a propagation constant, let us assume that P-j can be computed 
sufficiently accurately by the use of the following separable form: 


’ 1 = 1 


n=->oo 




-1(y, + 2rm) x/1 
e e 


(3.25) 
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Substituting Eq. (3,30) into Eq. (2.11) and simplifying, we obtain 

fj(z) + = 0 (3.26) 


where 

= C(a, - v2)(iil^/il)^ + + 2io)v(iy^/il)]/a^ (3.27) 

and ^n ~ ^ solution to Eq. (3.25) can be written as 

ik„z -■>k„z 

f„(z) = B,(n) e " + B^(n) e " (3.28) 

However, imposition of the radiation condition (i.e., P.j can propagate in 
the domain z < 0 only) results in 

B2(n) =0 

ik_z 

and f^(z) = B^(n) e " (3.29) 

In order to determine B-j(n), we use the B*C given by Eq. (2.12); namely 


3P^ 

3z 


= -D ( — + V — i W 

3t 3X ^ ” 


Z=0 


which gives 



= I i B^(n) exp {-iy^^ x/l} 

n=-oo 

ly x/A 

Multiplying both sides of this equation by e and integrating w.r.t.x 

from 0 to ly one obtains 


B-j(n) = w^. B^(n) 

where 
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V Pi 4 . (-0)^ + 2i0)VX. + V^X?) r A.£+iX 1 

' ’ V j=H (iu„/A+Xj ^ / 


(3.30) 


Combining Eqs. (3.30) and (3.35), the expression for the induced pressure 


P-j becomes 


ik:,z -ivA 


^ I B^(n) e " e " e 


(3.31) 


n=-oo 


and 




(3.32) 


Treating P-j as an external force to Eq. (3.2), the particular 


solution of Eq. (3.2) (call it Wp+) is 


W, 


p+ (x) = w^. Q+(x) 


where 


and 


-ip^x/A 

Q+(x) = I D/q+(n)] e " 

n=-oo 


q+(n) = [1/B^(n)] {D(y^A)'^ - A 


(3.33) 


(3.34) 


(3.35) 


Combining Eqs. (3,3) and (3.33), we obtain 


4 X.x 

& a ^ 


w+ < I a. e + Q.(x) ) *, +ve going wave 

1 j=i ^ / 


W(x) = < 


4 Xjx 


w < I a" e + Q (x) ) ; -ve going wave 

L-\d=i j ■ j 


(3.36) 
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I where Q_(x) has a similar form as that of Q^(x) with the .obvious changes 
which correspond to a negative going wave (i.e., replace y, at and w^ by 
-y, a^ and w_, respectively). 

Keeping in mind that the total solution should satisfy the boundary 
conditions at the extreme ends, the unknowns w^ and w_ are then determined 
from such conditions. Assuming that both ends are identical to the internal 
supports (i.e., transversely rigid, and have the same rotational stiffness 
k^), the total bending moment at both ends should be equal to zero; i.e., 

M^(0) = M^(HB) = 0 (3.37) 

where B is the number of the bays of the finite beam. The total moment at 
either end consists of the contributions of the forced response of the 
infinitely long beam and that of both the positive and negative going waves 
of the semi-infinite beams. In other words, 

M^(0) = 0 = M^(0) + M_(0) + (0) - k^ W|.(0) (3.38) 

where M.^.p(x) = contribution of the forced response to the bending 

moment at x. 

W^(x) = total slope at x. 

Note that the rotational spring k^ at x = 0 is not a part of the periodic 
element between x = 0 and x = il which explains the presence of the last 
term of Eq. (3.38). Similarly, 

M^(ilB) = 0 = M^(£B) + M_(AB) + M.^^(fi,B) (3.39) 

Substituting for the values of M^, M_ and in Eqs. (3.38) and (3.39), we 
obtain two simultaneous equations for w^ and w_ as to be explained below. 
Note that 
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4 . • 4 

W;(0) = w^.{ I a. X: + Q (0) }+ w { I aT X. + Q (0) } + W1 , 
^ j=l J J + j=T J J - ■ 1T>T 


( 0 ) 


and M(£B) = M(0). Thus, Eq. (3.38) becomes 


^ Q;«» - T < J, “J ^ 


k_ 4 


(0))} 


*W. a aj x2 ^ q:i( 0) - f a: X. + Q:(0))) = f w;„f(0) -u;„j(0) 


"r « 




and Eq. (3.39) reduces to 


w+ I A A* ^ I ® 

j-^1 J J + J J - 


= -W5„, (.B) 


Hence, from Eqs. (3.41) and 3.42), we find 


= 




®2“l ' ®1“2 


where 


6] w;.^^(0) - w;„f(0)} + w*(£b) 


w = 


^l“2 ■■ ^2“l 


k„ 4 


a, = 


«o = 


-I a; xf . q;(0) - ^ iX^ aj X. * q;(0)} "| 


I a: X? + Q"(0) - { I 

j=i ^ ^ ^ j=i ^ 

j=l ^ ^ 


kr 4 


3-j = e 
82 


4+0 iyyB 

■ e ^ q;(0) 






4 P -iyjB 
I a: XT + e Q"(0) 
j=l 


(3.40) 


(3.41) 


(3.42) 


(3.43) 


(3.44) 


(3.45) 
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Once Wj. and w_ are determined, the total displacement (which is the 
sum of the displacements given in Eqs. (3.36) and (2.18) with the term 
omitted) and other related quantities are readily determined. Since 
we have substituted a unit plane wave for the excitation pressure P, the 
computed response is actually the frequency response function H. Therefore, 
we can write 

H . X.x 4 X.x 

H(x,jp) = w^. { at e J + Q^(x)} + w_ { aT e ^ + Q_(x)} 


« -i(PQ+2nTr)x/£ -ipox/il 

I A {e . e } 

n=-oo 


0 < x < £ 


(3.46) 


Again, Pq is the forcing wave propagation constant and the coefficients 
A^'s are obtained from Eq. (2.33). 

It is to be noted that the frequency response function or any related 
quantities outside (0 £ x £ H) can be obtained most simply from Eq. (3.46) 
by using the wave property given in Eq. (3.5) (note that for a negative- 
going wave the term e“^ replaces e'^ in Eq. (3.5). 

B. Transversely Elastic Supports 

The analysis in this section is basically the same as that of the 
previous one (i.e., finite beam on rigid supports). However, there are 
some differences in the B.C.'s at the internal supports and in the number 
of free waves which may propagate in the structure. This will become clear 
as we proceed to calculate the coefficients A^'s of Eq. (3.3). 

Using the wave property as expressed in Eq. (3.3), we obtain for 
(i) the displacement 

W(il) = e^W(O) 


41 



or 


4 X.l „ 4 

I A, eJ .e'* J Aj 

j=l ^ j=l J 


(ii) the slope 

The same relationship as in Eq. (3.9). 

(iii) the moment 

The same relationship as in Eq. (3.11). 

(iv) the shearing force (referring to Fig. 5) 

Vi ■ ''t ° 


or 


W'"(Jl) + W(£) = e^W" (0) 


Substituting Eq. (3.48) into Eq. (3.3), 


I A,e«>‘= e*- I A a] 

j=l J J j = l ^ J 


(3.47) 


(3.48) 


(3.49) 


Combining Eqs. (3.9), (3.11), (3.47) and (3.49), we obtain a set of 

simultaneous equations relating the coefficients A-'s and the propagation 

%} 

constant y. This set of simultaneous equations can be written in a 
matrix form as 


G {A} = e^ Q {A} 


where {A}‘ = 


G = 


[A-j , A 2 , A^ , A^3 




X 3 A 

X.l 

e e ^ 

e 

e 

Xt ^ Xa^ 

X^e X^e ^ 

X^e^ 

X z 

V 


(3.50) 

(3.51) 


,2 . '^rx_^l^ 


O 

.C ^ C 


2 . '^rx.^3^ /,2 . '^r>_V 


(X^ + -f)e (x; + -ff)e ^ (X^ + -rf)e (X^ + -jf)e 


"2 D 


3 ^ D 


'4 D 


k. X,£ 

^ I 


kt X 2 JI 


(^T + -ir)e ' {X% + -rf)e (X 3 + -f)e (X^ + -jf)e 


k. X-Jl 




‘4 D 


(3.52) 
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station r station r+1 


Fig. 5. A Single Periodic Element 

(Transversely Elastic Supports) 
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1 


1 


1 




'3 

.2 





(3.53) 


= the spring constant of the elastic support, 

and T denotes the transpose of a matrix. A matrix of the form given in 
Eq. (3.53) is called a Vandermonde matrix. Since the roots X.'s can be 
assumed distinct, it can be shown that this matrix is nonsingular and 
therefore exists [23]. 

The remainder of the analysis can be carried out parallel to that of 
the previous section. Rewrite Eq. (3.50) in the form of 


IT{A} = e^{A} 

where 

C = q"^ G 


(3.54) 

(3.55) 


The matrix C is cross symmetric and possesses two pairs of reciprocal 
eigenvalues [16]. If these eigenvalues are denoted by A.(j = 1^4), 

J 

and let 



(3.56) 


then the corresponding two pairs of propagation constants are 

Pi = logg = -P3 

Vg = log^ Ag = (3-57) 
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In other words, there will be two free waves propagating in the positive 
direction with the propagation constants and P2» and two free waves 
propagating in the opposite direction with -y-j and -pg. Also, the 
eigenvectors corresponding to the different propagation constants can be 
determined. Of course, only the normalized eigenvectors (the components 
of an eigenvector divided by one of them) are unique. If we denote A./A, 
(corresponding to pj^) as |^ , then the displacement Wj^ due to the free 
wave with propagation constant pj^ can be written as 

4 X -X 

Wk(x) = Wk I a^^k ® ^ (3.58) 

j 1 

where Wk is the amplitude of the wave and will be determined from the 
B.C.'s at the extreme ends of the finite beam. 

The effect of the induced pressure P-j can be accounted for in a 
manner similar to that explained in the previous section. To do that, 
let P-j k tie the pressure induced by the wave with propagation constant Pk , 
and let it be expressed as follov/s: 

Pl,kl ^ \ exp{-(pj^l^ + 2mr) y./l] e''“^ (3.59) 

2=0 

where B-j k(f^) given by Eq. (3.30) and Pj k the imaginary part of Pk* 
Similarly, let the particular solution of Eq. (3.1), with P-jj ^ replaced 

p ^ 

by P-J k| » be Wk* It is simple to show that 

* p 

Wk = Wk ITk (3.60) 

CO ^ ”''^^I,k 

where W. = ^ B, . (n) — ^ (3.61) 

K n=-» p, . + 2mr 4 , 

) - Cg J 
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We shall now proceed to determine the unknowns using the B.C.'s 

at the two extrane ends of the entire beam. Let the supports at the 
extreme ends be Infinitely rigid In the transverse direction, and let their 
rotational stiffness be k^. Expressed mathematically, these conditions are 


4 p 

Wt( 0 ) = I {w|;( 0 ) + W^( 0 )} + W.^^(O) = 0 ( 3 . 62 ) 

Ic*”! 

W^(£B) = I {w]^(ilB) + \{SLB)} + W.^^(aB) = 0 ( 3 . 63 ) 

k”” *l 

4 

\{0) = I fT^(O) + - k^ W-( 0 ) = 0 ( 3 . 64 ) 

k“l 

4 

M^(AB) = I M|^(AB) + = 0 ( 3 . 65 ) 

k*™"! 


p 

where TT|j(x) Is the moment corresponding to {Wj^(x) + Wj^(x)}, and the 
subscript t denotes the total value (I.e. , total displacement, moment, etc.). 
Application of Eqs. (3.62)-(3.65) results In a set of simultaneous equations 
for the unknowns Wj^'s. These equations are combined Into a matrix equation 
as follows: 


where 


V{w,^} = {R> 


{R}^ = [-W.„f(0), -W^pfUB),-^W!^f(0) 


(3.66) 


WV^f(O), -W5^^(£B)3, (3.67) 


and the elements of the matrix V are given by 

4 

,(0) 


'l,k = ®j,k ^k< 


^2,k = ® 


Pi-B 4 ^ ipj i^B 

" I a + e a (0) 

j = l 
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''3.k 

''4.k 


jl, ®j.k “k^^^ " D ^ ®j,k 

4 yv o iy-r |-B 

® jJ, ?j.k 4 ^ *'k"’> " 


"r ^ 


k = 1 ^ 4 


Once the unknowns W|^'s are determined, the total displacement is 
readily determined by summing up the displacements corresponding to the 
free waves W|^*s and that of the infinite beam. Since we have substituted 
a unit plane wave for the excitation pressure P, the total displacement is 
actually the frequency response function H. Therefore, 


‘ {Wk(x. 


(3.68) 


where Hinf(^* 9lven by Eq. (2.60). 

Again, the frequency response function or other related quantities 
outside (0 _< X £ £) can be easily' obtained from Eq. (3.68) by using the 
wave property given in Eq. (3.5) (note that for a negative-going wave, 
e"^ replaces e'^ in Eq. (3,5) ). 
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IV. TWO-DIMENSIONAL CASE: AN INFINITELY 
LONG PANEL BACKED BY AN AIR CAVITY 


In the preceding sections, some possible designs of passive compliant 
surfaces have been analyzed using simplified one-dimensional models of 
periodic beams. Such simplified analyses reveal basic characteristics of 
real structures, but more accurate representation requires the use of two- 
dimensional models. 

In this section, an infinitely long panel backed by an air cavity of 
depth d is considered. The side walls of the cavity are assumed to be 
either acoustically hard or acoustically soft and both cases will be dealt 
with in this investigation. The panel is supported by pre-tensioned 
circular solid wires at distant I apart. The tension in each wire is T. 
The configuration of such a panel is shown in Fig, 6. 

The equation of motion of such a panel, not directly over a wire, 
is given by 


where 


DV^W + mpW = 




^ 2 ^ I 


( 4 . 1 ) 


v4 2 ... s'* + 

'' 4 ^ 22 4 

8x 3x 3y By 

P = the turbulence pressure field (if the panel were motionless) 
P,| = induced additional pressure generated above the panel 
?2 = pressure variation transmitted into the cavity 
mp = mass per unit area of the panel 
D = flexural rigidity of the panel. 


Again, we shall determine first the frequency response function. For this 
purpose, P is replaced by the sth plane wave gisiry/b^ 
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Fig, 6 An Infinite Panel Backed by an Air Cavity 


The panel is assumed to be simply supported at both .sides y = 0,b. 


Accordingly, the boundary conditions for W are 

W(x,0,t) = W(x,b,t) = 0 (4.2) 

The pressure which is induced by the panel motion is governed by 
the convective wave motion 

" '' '■i ■ 


This induced pressure, , is subjected to the B.C. 


a a ? 

^1 ^ 3x^ “ 3z 


1 


z=-d 


(4.4) 


and the radiation condition that it propagates in the domain z < 0. 

The pressure wave ?2 which is also generated by the panel motion and 
transmitted into the cavity is governed by the equation 



2 2 


(4.5) 


where ag is the speed of sound in the cavity. If the side walls of the 
cavity are acoustically hard, then the boundary conditions for are 


aP 

dz 

3P 



at z = 0 
along y = 0,b 


3P, 


3z 


z=-d 



(4.6) 

(4.7) 

(4.8) 


where P 2 is the air density inside the cavity. 
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Since the displacement must be spatially periodic (i.n the x-direction) 
and must satisfy the B.C, (4.2), it is convenient to write the displacement 
W due to the excitation e^^“^ “ gisTry/b follows: 


«,(x.y.t)= I, I A^„ e-'V/*' sin (SSy) 


m=l n=-~ 


(4.9) 


where ^0 = ^ 0 ^ ^0 forcing wave propagation constant (= k£). 

Again, our immediate objective is to determine the coefficients A „'s. 

mn 

In view of the B.C. given in Eq. (4.4), the induced pressure P-j may 
be expressed in a form similar to that of the displacement; namely 


P,(x.y.z,t) = I e-’V/* sin (!^y) e^“‘ 


m=l n=-«> 


(4.10) 


Substituting Eq. (4.10) into Eq. (4.3) and simplifying, one obtains 

f" ( 2 ) + kL (z) = 0 

mn mn mn ' ' 


(4.11) 


where 


af a^ £ a^ £ £ b 


(4.12) 


The solution to Eq^ (4.11) can be written as 

f (z) = Bi(m,n) e’*Snn^ + 82 ( 01 , 0 ) e"^*^mn^ 


(4.13) 


Since P-j can propagate in the negative z-domain only, B 2 (m,n) must be equal 
to zero and Eq. (4.13) reduces to 


ik_z 


f^^(z) = B-|(m,n) e mn 


(4.14) 


Using the B.C. given by Eq. (4.4), 
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B^(m,n) = 


-ip^u)a^(C^ - V)^ 

r tin \i\^ / BITT i. >2i^2 

C„{(C„ - V) - a, - (-5- • ^a,) )“ 


mn 


(4.15) 


where C = ^ . If we substitute Eqs. (4.14) and (4.15) into Eq. (4.10), 


n y, 


we obtain 


P^(x,y,z,t) = I I B,(ni,n) e’*^mn^^‘''^^ e"’^n^^^ sin (^^y) e^“^ (4.16) 
m=l n=-» 


and 


P,(x.y.z.t)l . I I e-’Pn’</'‘ sin (4.17) 

2 =-d m=l n=-a> 


where B|(m,n) = B^(m,n) e“'^mn°^ It is to be noted that the comments which 
have been made concerning the effect of P.j on the response are the same for 
both the one-dimensional and two-dimensional cases (see page 16). 

Similarly, because of the B.C. (4.8), a convenient form for ?2 is 


Pp(x,y,z,t) = I I 


m=0 n=-«> 


mn 






cos (¥-y)e 


j‘(Dt 


(4.18) 


The B.C.'s given by Eq. (4.7) are automatically satisfied. Substituting 
Eq. (4.18) into Eq. (4.5) and simplifying, one obtains 

^ ’in = 0 

“hore kL = 4 - (X)' - (t>' • <*•“> 

The solution to Eq. (4.19) may be written as 

Fjjjj,(z) = D-|(m,n) e^*^mn^ + D 2 (m,n) e'^’^mn^ (4.21) 
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Imposition of Eq. (4.6) results in 

D-j(m,n) = D 2 (m,n) h D(m,n) 

and ^mn^^^ 2D(m,n) cos (4.22) 


The B.C. given by Eq. (4.9) leads to 


1 1 ^mn cos (^y) D(m,n) K sin (K d) 

„ mn ' b mn ' mn 


m=0 n=-“ 




(4.23) 


Multiplying both sides of Eq. (4.23) by e’^n'^/^cos (^^^ y) and integrating 
over (0 < X < A; 0 < y < b) , 

D(m',n’) K^.„, sin (K^,„, d) 


= Pg 03^ A I [ Sin (2f y) cos (S^|^ y) dy 

m=l •'r, 


which becomes, after simplification. 


Vn' ' 


2 P2 u) 


L ^mn' 


2m 


f(m') K^.„. sin (K^,^, d) m=l - m’^) 


(4.24) 


where f(m')=^ 


11 ; m’ 0 

12 ; m' = 0 . 


(4.25) 


Substituting Eq, (4.24) into Eq. (4.18) gives 


m=0 n= 


2 ~ » 4 cos (K_ z) 


mn 


•iPpX/A 


fO n.- f(m) . sin (K,„„ d) 
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(4.26) 


cos 


I 

r=l 

m+r=odd 


'itot 


rn 


ir(r^ - 


m^) 


Exchanging the symbols r and m, and reversing the summation order in 
Eq. (4.26) to make it compatible with the form of the displacement W, 
the expression for becomes 


P 2 (x,y,z,t)l P 2 I J e \ 

V r+m= 


4m 


=0 f(r)iT(m^-r^) 

Fodd 


cot (K^^d) cos (^y)/Krnf ® 


imt 


(4.27) 


It is to oe noted that when is imaginary, the quantity 


cot (K^^d) 


becomes - 


Coth (iK^Jd) 


rn 


IK, 


. This means that the nth component of Pg would 


rn* 


act as an added mass. 


Determination of the coefficients Aj„p's: 

The virtual work principle will be used to determine A „'s. The 
total virtual work done by all the forces on and in the panel (including 
those due to the wires) should be equal to zero. To calculate the virtual 
work, we will consider one periodic element only for the same reasoning 
given for the one-dimensional case. Let the virtual displacement be 

sin (^ y) (4.28) 


As usual in complex algebra, the complex conjugate of 6W will be used in 
calculating the virtual work. 
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Let the virtual work done 

rbf£ «> » 


by + mpVi - I J be 


5W 


1 


rDfJO » «> 

LL i ^ 

JO'to m=l n=- 


e-'fV'n')’'" A 

tn n 


mn 


sin y) sin y) ^(m.n) dx dy 


»m'Tr 


where 


4>(m,n) 


,n) - {(x>* + 2( 

i p, <«a,(C„ - V) 


m/tr . ^n\2 . #mir\4 


f + } D - mpo,^ 


b I' ■ '~b 
2 


) 


- «1 - (^1 X 

Carrying out the integration in Eq. (4.29) and simplifying, we 

«'fl =4Vr' «Vn' 

Similarly, the virtual work done by is 

= [X 


2=-d 


= 6A„. , A y A„„, 
m n mn 

m= I 


P 2 <ii m 4 cot (Kppd) 


-5 5~ ' 

Tr(m^ - r^) f(r) K, 


rn 


j'’ sin (S!> y) cos (if y) dy 


I •■J-, ''mn' “m'n' 

m=i 


where 


(4.29) 


(4.30) 

obtain 

(4.31) 


(4.32) 
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ii»(m,m'.n) = I 


y 

4 pg w tran' 4 cot 


r=0 (m^-r^)f(r)Ti^ 'K, 


m+r=odd 

m'+r=odd 


rn 


(4.33) 


The virtual work done by the excitation 


6Wp = 


A„,„. sin (2^y) 


gisiry/b ^jy = (s,in',n') 


where 


(4.34) 


-I (m'TT/b) {1 - 




2 


a{s,m',n‘) = < 


(V) 


(f) 


^ *^\'n' 2 “^n'O 


6^.0 sfm' 


; s = m' 


(4.35) 


The wire, because of its elastic and inertial properties, will produce 
shearing forces and moments in the panel along the line of attachment. The 
virtual work done by the shearing forces is (say at x = 0) 


^ fb 


^W.s = 


8^W 


9 2 

9 W T 9^Wi 


^w T ^ pw ''w ^ dy 

W W W W 


-b OO oo f 
= ^ 

•'o m=l n=-» ^ 


I (li^) _ p A u + T( 

w ' b ^ ^w w 




\n “m'n’ <T 


/m’lT 


= 1 ^ %-n ^^m-n- 


n=-o 


(4.36) 


where 
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Ey = Young's modulus of the wire materials 

I., = moment of inertia of the wire about the x-axis 
w 

= density of the wire material 
= cross sectional area of the wire 
T = tension applied to the wire 
and 

C(n.',t.) = (4.37) 


Also, the virtual work done by the torsional moments in the wire is 




’w ■>>. W 0 


6W'(0,y) dy 


(4.38) 


where G is the shear modulus, C is the Saint-Venant constant of uniform 

torsion, and is the polar moment of inertia of the wire. Carrying out 
w 

the integration in Eq. (4.38) and simplifying: 



I 


n=-oo 


^m'n 


6A. 


m'n' 


r(m' ,n,n') 


(4.39) 


where 

r{m'.n,n') = (-GC{iS^)^ + (4.40) 


Since the total virtual work must be zero. 


6W, + 6W_ + 6W^ + 6W„ + 5W^ = 0 
I P» s m p 


(4.41) 


m 


Substituting into this equation the values of 6¥-|, 5Wp , 6Wp, and 6W^ 
from Equations (4.29), (4.32), (4.34), (4.36) and (4.39), respectively, we 
obtain a set of simultaneous equations for the coefficients 
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and 


Pptx.y.z.t)! = I I cot (ic_d) 

z=-d m=l n=-<» ^10 


nw „\ ^iwt 




(4.54) 


Now, the virtual work done by the pressure ?2 ’’s 

rbr£ ~ “ P-, (u2 


6l 


P2 




h P7 

* ^ 2 * k;;; 7;;7 \'n' 


(4.55) 


Therefore, using Eq. (4.55) instead of Eq. (4.32), we obtain a new set of 


simultaneous linear equations for the coefficients 

- 1 " - 
A^.„i 4>(m' ,n‘) + -^ I A(m' ,n,n') «e(s,m',n') 


n=-oo 


where (|)(m’,n') = <|)(m',n') + P 2 u)' 


2 

•Sii'n' 


(4.56) 

(4.57) 


and 6(s,m',n') is given by Eq. (4.43). Moreover, we can write a truncated 
Eq. (4.56) in a matrix form 

P^nc.6^ 'W ' <B{s,m,n)} (4.58) 


where 


^naB = (j)(m,n) + ]; A(m,3,n) , 


(4.59) 


m = 1 ,2 N. 

0=1 ,2 , . . . , N^ 
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n “ j **1 (0(1 j 

0 “ “N^( • • • ( ”1(0(1( . . ( ( 

and (2NJ +1) and (2N^ +1) are the number of significant modes in the y 
and X directionS( respectively. 

Once the coefficients are determined from Eq. (4.42) or Eq. (4.57), 

the frequency response function of the structure due to the sth component 
of the pressure wave, is readily determined; namely 


H (x.y,k) ^ I I 
* m=l n=- 


mn 


-ip x/A _ 

® " sin(!^y) 


(4.60) 
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V. STRUCTURAL RESPONSE 

In the preceding sections, we were mainly concerned with the frequency 
response function, H. The next question is whether or not desired structural 
motion can be generated by boundary layer turbulence excitations. Under 
such excitations, the structural response is a random process and should be 
characterized in terms of statistical quantities such as spectral density, 
mean square value, etc. 

-If we accept the frozen pattern hypothesis for the turbulence field, 
then after reaching stochastic stationarity the displacement W may be 
expressed as 

I H5(x,y,k) fg e’“^ dF(k) (5.1) 

—00 S=-00 

where E{dF(k^) dF*(k^)} ^ S(kj - kg) Sp(k^) dk^ dkg 



and f is given by Eq. (1.4). The frequency response function H is 
given in Eq. (4.60). The autocorrelation function of the displacement W 
is 

r„(x,y,t^,tg) = E{W(x.y,t^) W*(x,y,tg)} (5.2) 

'V 

Combining Eqs. (5.1), (5.2) and noting that the Fourier components f^ are 
uncorrelated, we get 

R„(x.y,T) = r^(x,y,t-,-tg) = T J |Hg(x,y,k|^ f^ 

J CSwOO 
••00 ^ 

Sp(k) e^'^^V dk (5.3) 
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The spectral density of the displacement, is the Fourier transform 
of the autocorrelation function; i.e., 

*„(x.y.J2) = ^ I" Rv,(x.y,T) dx 

«aOO 

««00 —00 
gi(kUj. - 

= IT r ^ lHgCx.y.jf^)l^ f S (jf-) 5(fi -w) doi 
‘'c S=-“ c “c 



(5.4) 


For the one-dimensional case, becomes 

° ^p<“> 

0 


(5.5) 


If the drag reduction mechanism hypothesized in Reference [31 is 
true, then a good design is one that can respond to the boundary layer 
excitations in a desirable manner within a particular frequency range. 
The root mean square, value of the structural response attributed to a 
frequency range is ; 


W 

r.m.s 


r* W,i 


(x,y) = 




^„(x,y,oj) d(u 




(5.6) 
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where ui^ and are the lower and upper limits of the frequency range, 
respectively. 

To carry out the integration in Eq. (5.6), we need an expression for 
the turbulence spectrum. Considerable amount of information is available 
on experimentally determined boundary-layer turbulence spectra at various 
free stream velocities. See, for example, [21,24]. In particular. Bull 
[24] has done extensive spectral density measurements for subsonic boundary 
layer pressure fluctuations. Therefore, numerical calculations in this 
thesis are based on Bull's spectrum but with the decay factor in the 
convection direction omitted, and the results will be presented in the 
section dealing with the numerical results. 
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VI. SPECTRAL DENSITIES OF AND (P^ + Pj 
A. Spectral Density of P-j 


Let the frequency response function of P-j due to the excitation 
gisiry/b gi(o)t - kx) _(x,y,k). After reaching stochastic stationarity, 

I 9S 

the pressure P-j can be expressed as 

Pl(x,y,t) = r i e^“* dF(k) (6.1) 

««-00 ^ 

Carrying out an analysis similar to that of the displacement, the spectral 
density of P^ (for the two-dimensional case) can be expressed as 



(x,y,ai) 



( 6 . 2 ) 


where ^(x, jp) is obtained from Eq. (4.17), namely 
0 


00 00 ^ • i ft 

Hi s(x.y, if-) * I I B,(m,n) sin (^ y) 

“c m=l n=-«> ' ^ 


(6.3) 


For the one-dimensional case, 

$p^p^(x,a)) = lH^(x, ^)|^ Sp(co) (6.4) 

In this case H-j depends on the configuration of the beam. For example, 
in the case of an infinitely long beam on rigid supports, H.j is given by 
Eq. (2.19). 

B. Spectral Density of (P + P]) 


Let 

P = P + P-J 


(6.5) 
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The autocorrelation function, »t 2 ) , of P is given by 

TppCx.y.ti.tg) = E{P(x,y,t^) PCx.y.tg) } (6.6) 

= E{P(x,y,t^) PCx.y.tg)} + E{P^(x,y,t-|) Pi(x,y,t 2 )} 

+ E{P(x,y,t^) P^(x,y,t 2 )> + E{P(x.y,t 2 ) P-,(x,y.t,)} (6.7) 

When the stochastic stationarity of P and P^ is reached, Eq. (6.7) becomes 
R^p(x,y.T) = Rpp(x,y,r) + Rp^p^(x,y,r) + Rpp^(x,y,x) 

+ Rp^p(x.y.T) (6.8) 


The spectral density of p is merely the Fourier transform of 
R^^(x,y,T); i.e., 

*pp{x,y,m) = $pp(x,y,oj) + $pp^(x,y,w) + 2Re{$pp^(x,y,o))} (6.9) 

Expressing p and p-j as in Eqs. (1.5) and (6.1), respectively. 



Rp^p(x,y,T) = E < J J I J _ e’'=2>‘ 


i(klti - kptn)^ ic 

e ‘ ^ dF(k^) dF (k 2 ) dk^ dkg 


( 6 . 10 ) 


‘ f J Sp(k) dk 

-00 S--“ 

(6.11) 

Taking the Fourier transform of Eq. (6.11) and simplifying. 


(x 

PPl 


,y,o)) = e’ u“ fg Sp(w) (6.12) 
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pp'^»y,w; = r - 

S=-co s 




(6.72) 


%(^oy,io) = 


f=’' 


3nd (6 ]■}) 

‘°-'3). we ffnt, 


i=>« ^ 


that 


<0 


2Re I e e^'isiry/b 


^or the 


one- 


^^mens 




Tonal 


one obtafns 

V'--) = S f„, f, 

«- +22efptT' 




(6.73) 


(6.14) 


(6.15) 
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VII. NUMERICAL RESULTS 


In the preceding sections, we have derived expressions for the frequency 
response functions and the root mean square values of the displacements for 
a number of structural configurations. In this section, a parametric study 
will be presented on the effect of the changes in span length, plate thick- 
ness,..., etc., on the structural response as represented by the change in 
frequency response function and the root mean square value. In order to 
carry out the computations, the pressure spectrum Sp(o)) has to be specified 
for which we chose Bull's spectrum as shown in Fig. 7. The coordinates of 

p 

the figure are Sp{a)) V/qQ 6* and a>6*/V where 

V = free stream velocity 

1 2 

qQ = dynamic pressure = ^ pV 

6* = displacement thickness of the boundary layer 

The displacement thickness is related to the boundary layer thickness 6 by 

6* = 6/8 (7.1) 

and according to Schlichting [25] 

6(x) = 0.37 X R“^/^ (7.2) 

where x = distance from the leading edge to the location 

where 6 is to be evaluated 
R = Vx/v^ = Reynold's number 

a 

V, = kinematic viscosity of the fluid medium (air). 

a 

Traditionally the boundary-layer thickness is measured to where the velocity 
is equal to 0.99 V. 


68 




Bull has fitted his experimental data by the expression 

5* — _ — 

Sp(u)) = {3.7 e'^“ + 0.8 3.4 }x 10"® (7.3) 

where u = aj6*/V (7.4) 

Thus, for a certain flow velocity V, 6* can be determined and in turn Sp(u) 
can be computed. 

Since the vibration is random, a measure of its magnitude is the root 
mean square value which takes the place of the amplitude of a deterministic 
vibration. Also, there is lack of a universally agreed upon notion 
replacing the usual concept of wave length. However, we shall propose 
the term dominant wave length as a measure of the spatial variation of 
the total vibration which is computed from the dominant peak frequency 
of the response spectrum. 

The numerical results presented hereafter have been obtained by using 
a double-precision procedure on the IBM 360/75 digital computer. Some 
brief comments about the computer programs involved are discussed in 
Appendix I. These numerical results cover: (i) the one-dimensional and 
(ii) the two-dimensional cases. 

(i) The One-Dimensional Case 

The structure which is idealized as a periodic beam is assumed to be 
made of mylar with the following properties: 
p (density) = 1.3915 x 10^ k /m^ 

D = Eh^/12(l -v^) 

V = .3 (Poisson's ratio) 
h = beam thickness (to be determined). 

E = 2.757 X 10^ N/m^ 

“ Ph (mass per unit area) 
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The physical constants of the air above the beam are 

Pi • 0.1227 k /m^ 

-8 

a-| ■ 340.036 m/sec 

V, ■ 2.88 X 10"^ m^/sec 

a 

A parametric. study using the above values has been carried out to 
obta1ii/.dh optlM^^^ In terms of the following physical quantities, 

which meets the requirements set forth In reference [5]: 

(a) the periodic length I 

(b) the beam thickness h 

(c) the rotational stiffness 

(d) the foutijdt^eri material properties (k^. k 2 > n^) 

(e) the elasticity of the supports k^ (In the case of elastic support) 

As explained In the text, the one-dimensional case Includes both 
Infinite and tilbtte. bsAms having supports which can be assumed transversely 
rigid or elastic. In all of these cases, the computations are made for the 
mid-span point (x ■ t/2). Also, the convection velocity U^. is taken as .8 
of the flow velocity V which Is chosen as 30.48 m/sea. The different 
numerical results, corresponding to different cases, are presented in the 
following subsections. 

A. An Infinite Beam on Rigid Supports 

Our parametric study has been restricted within the following ranges 
of values: 

t » ( 2.54 n, 6.35) x lO'^ m 
h - (5.08 'V 25.4) X 10"^ m 
k^ = (2.714 n. 13.57) X 10^ N/m^ 
kg * (2.714 X 10^ n, 2.714)x 10^^ N/m'^ 
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= 0 2.714 X 10^ N-sec/m^ 

k ■ 4.448 X 10 ^ 4.448 x 10 ^ N-m/m/rad. 

r 

The results of our study are described below: 

(1) A typical response of the beam, represented by the square of the 
absolute value of the frequency response function, is shown in Fig. 8. It 
can be seen from this figure that there is a dominant peak occurring at a 
certain frequency, say f„,„. The peak region is the dominant part of the 
structural response and, therefore, our attention will be focused on it. 

For the relatively low flow velocity (V = 50.48 m/soc) considered 
herein, the peak response frequency appears to be near the fundamental 
natural frequency of the system. An estimate of this natural frequency 
can be obtained by considering a simply supported beam of length I resting 
on an elastic foundation of stiffness (i.e. , kg*®). The natural 
frequency of the nth mode is given by 

f . 1 ^ (7.5 

n ZiT ' m^j 

If we use the values: 

jl = 3.81 X 10*'^ m 
h = 5.08 X 10'^ III 
k = 2.714 X 10^ N/m'^ 

Kl 

then the fundamental natural frequency f^ is approximately equal to 990 Hz. 
The peak frequency f computed by using the above t, h and k, values as 
seen in Figs. 8 and 10 is 900 Hz. The slight discrepancy between 990 and 
900 is due to the fact that other factors which affect f ^ have not been 
included in Eq. (7.5). Such factors are the convection velocity, the 
rotational stiffness, and the interaction between neighboring spans. 
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(2) Within the range of parameters investigated, the dominant wave 
component is found to be associated with an n value between -2 and 2. 

Some samples of the coefficients A^'s indicative of relative importance of 
different wave components are shown in Tables 1-3. It is important to 
mention that these A^^'s were computed from Eq. (2.33) or its simplified 
versions and by truncating the number of simultaneous equations to 16 
(i.e., N *= 8), It was found from the computations that a minimum of N = 8 
was necessary to obtain a reasonable accuracy especially at or near the 
peak frequency. 

(3) The change in the periodic length, Z, greatly affects the 
structural response including its amplitude, peak frequency, and the 
dominant wave length. As expected, an increase in Z generally leads to 
a reduction in the peak frequency and an increase in the amplitude. The 
wave length of the nth wave component is related to the periodic length Z 
as follows: 



n 


+ 



max 



(7.6) 


where we assumed that this component contributed to the response at the 
peak frequency f . Since the dominant terms are those corresponding to 
n between -2 and 2, the design criterion concerning the wave length should 
be checked against the wave length of each of these components. For 
instance, when n = +2, the absolute value of the R.H.S. of Eq. (7.6) is 
approximately equal to 1/2 (since the value of f Z/\} is small within 
the range of values being considered). Therefore, the corresponding wave 
lengths 2 approximately equal to ZfZ. Consequently, these two 
wave lengths can easily satisfy the wave length requirement (£ 2.54 x 10 m 
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for values of I ranging from (2.54 to 3.8I) x tik> term wiiich poses n 

great difficulty in meeting the wave length requirement is A -j. The wave 
length associated with this term is obtained from Eq, (7.6) (by putting n 
equal to -1): 



The absolute value of the R.H.S. of Eq. (7.7) is always greater than one. 

Therefore, to obtain the desired wave length, a has to be less than 2.54 

X However, such a small £. would increase f beyond the desired 

max 

frequency range. As a compromise we have chosen Z as ^ lu \ for this 

parametric study (which might not satisfy the wave length requirement but 
it keeps f^^^^ within the desired frequency range). 

(4) The thickness h has a great effect on the amplitude and on the 

peak frequency of the response. An increase in h shifts to a higher 

value and generally leads to a decrease in the amplitude. 

(5) The increase in the stiffness k^ reduces the amplitude and 

increases the peak frequency (See Fig. 9) 

(6) The change in the stiffness V .2 does not affect the peak frequency 
appreciably. However, an increase in k 2 increases the amplitude. The 
limiting case of an infinite k 2 is simply that of an elastic foundation of 
stiffness k-j. The effect of the change of k^ on the response is shown in 
Fig. 10. 

(7) The increase in the damping coefficient for the range of 
values used, leads to an increase in the amplitude without any significant 
change in the peak frequency. (See Fig. 11) The structural response 
corresponding to the special case of = 0 is shown in Fig. 12. 
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(8) An increase in the rotational stiffness increases the peak 
frequency and slightly increases the amplitude. (See Fig. 13) 

Some representative root mean square values of the response, and the 

wave lengths of the dominant components are shown in Table 4. As lean be 

! 

seen from this table, the root mean square value, W„ „ ^ , is smaller than 

r»ui» s • i 

the desired values for most of the cases shown. Also, the wave length X_^ 
is always higher than the desired value for the reasons we explained before. 

It should be noted that, unless stated otherwise on the figure itself, 
the following values have been used in computing the response shown in 
Figures 8 through 13: 

S, = 3.81 X 10“^ m 

; 

h = 5.08 X 10 ^ m I 

= 2.714 X 10^ N/m^ 

k^ = 2.714 X 10^ N/m^ 

kp = 4.448 X 10"'^ N-m/m/rad. 

= 2.714 X 10^ N-sec/m^ 

= .01 

Hb = 0 

B. An Infinite Beam on Transversely Elastic Supports 

The only difference between the structural configuration of the beam 
considered here and that of an infinite beam on transversely rigid supports 
is in the introduction of transversely elastic springs, each characterized 
by a spring constant k^, at the supports. Therefore, in this section, we 
are mainly concerned with the effect of these transversely elastic supports 
on the structural response. Such effect is shown in Figs. 14-17. Figure 14 
shows a typical response of the beam over a wide range of frequencies, 
while Figs. 15-17 focus on the peak region only. Based on these figures. 
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we conclude that: 

(1) The elastic spring constant, k^, has a considerable effect on the 
location of the peak frequency. As k^ increases, the peak frequency is 
shifted to a higher value until it approaches that of the transversely rigid- 
support case (i.e., k^ = “). Thus, by choosing a small k^, it is possible 

-3 

to reduce the periodic length to as small as 2. 54 x lO m so that the require- 
ment for wave length can be satisfied without raising to outside the 
desired frequency range. 

(2) The amplitude of the response increases as the value of k^ 
decreases. Therefore, a suitable choice of k^ can help in bringing the 
amplitude to the desired value. 

Typical values of the root mean square displacements and the wave 
length of the dominant terms are shown in Table 5. It can be clearly seen 
that, except for the first one, the amplitude and the frequency of all cases 
presented fit the desired ranges. Moreover, the wave length X_-j of the 

_3 

component n = -1 is very close to the desired value of 2.5 x lO m (the 
difference is approximately 1%). 

It is to be noted that the following values have been used in computing 
the response shown in Figs. 14-17: 

£ = 2.54 X 10 

h = 5.08 X 10"^ra 

k = 4.448 X lO"^ N-m/m/rad. 

r 

n = 2.714 X JO^ N-sec/m^ 

'v 

= .01 
= 0.0 
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Table 1. Sample of the Coefficients is (in meters) 


n 

' -8 

-7 

-6 -5 

-4 -3 -2 

-1 ' 

0 * 

Re{A .} 
n 

4.26 X 10 

“^6.35 X 10“^ 

il.03 X 10“'^i.85 X 10“^ 

: 

3.88 X 10“"^!. 19 X 10“^ -1.06 x 10^ 

-1.75 X 10“^ 

1.29 X 10 

Im{A } 
n 

|-2. 26x10“ 

^ -3.43x10“^ 

• 

-5.58x10“^ ;^1. 03x10“^ 

-2.26x10“^ -7.5x10“^ 5.5x10“^ 

- ; 

1.5x10“^ 

-8.09x10“ 

— 

n 

J 

1 - 

i 

; 2 

! 

3 ; 4 

5 6 7 

8 

J 

• 

— — 

Re{A } 
n 

r 

i _•? 

16.6x10 ^ 

1 

■ 1 ■ 

1-7.44x10“^ 

' 

-9.65x10“'^ 1 -3.63x10“"^ 

-1.78x10“^ -1.01x10“'^ :-6. 09x10“^ 

i 

-4.19x10“^ 1 


Im{A } 
n 

1 -1.75x10“ 

{ 

^ 3.11x10“^ 

: 

4.04x10“^ 1.54x10“^ 

} . 

7.77x10“^ 4.47x10“^ • 2.84x10“^ 

i 

1.9 X 10“^ 

- 



i = 3.81 X 10 

mj h = 5.08 X 10 ^ m; 

= 2.714 X 10^ N sec/m^ 




= 2,714 X 10^ N/ra^; = 2.714 x 10^ N/m^; k^ = 4.448 x lO"^ N-m/m/rad; n^. = .01 
f = 898 HZ (peak frequency) ; * dominant term 


Table 2. Sample of the Coefficients (in meters) 
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Table 3. Sample of the Coefficients A^'s (in meters) 
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Table 4. Sample of Displacements and Wave Length of an Infinite Beam on Transversely Rigid Supports 
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Table 5. Sample of Displacements and Wave Lengths of an Infinite Beam on Transversely Elastic Supports 
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C. Finite Periodic Beam 

In this section, numerical results of the structural response of a 
5-span periodic beam will be presented. These results correspond to two 

T' 

cases: (a) the interior supports are transversely rigid, (b) the interior 

supports are transversely elastic. ! In both cases, each of the end supports 
is assumed to be rigid in the transverse direction, and to possess a 
rotational stiffness k^. The physical data used in the computations for 
these cases are the same as those used previously for the corresponding 
infinitely-long beams. Also, the* response referred to is for the midpoint 
of the first span. 

(a) Transversely Rigid Interior Supports 

Comparisons between the structural response of a finite beam and that 
of an infinitely long beam, other physical data being identical., are 
presented in Figs. 18-20. Each of these figures shows the square of the 
absolute value of the frequency response function vs. the frequency. 

Figure 18 shows the comparison when a viscoelastic foundation ik used for 
both finite and infinite beams, while in the case of Figs. 18 ahd 20 the 
foundation is purely elastic. 

As expected, the response curve for each finite beam contains multiple 
peaks clustered in frequency zones. Within the frequency range' shown in 
the figures, the response peaks are clustered roughly in the first free- 
wave propagation zone [15]. The number of peaks in each zone is expected 
to be equal to the number of spans of the finite periodic beam.- However, 
it is possible that some of those peaks might not appear because of the 
amount of damping present in the system or because of the choice of a 
coarse frequency interval in the computations. ^ 
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Fig. 20 Response of a Periodic Beam (Transversely Rigid Supports, 0) 
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(b) Transversely Elastic Interior Supports 

In this case, there are four free waves which might propagate in the 
finite beam, compared to only two free waves in the case of transversely 
rigid-support case. This increase in the number of waves, and the ill- 
conditioned matrices involved make the computations more complicated. More 
specifically, the Vandermonde matrix given in Eq. (3.53) and the system of 
simultaneous equations relating the unknowns Wj^'s (i.e., Eq. (3.66)) are 
very ill-conditioned (meaning a small error in the elements of these matrices 
would magnify leading to much larger errors). Although there are special 
techniques to deal with such matrices (see Appendix I), the computational 
accuracy can suffer greatly from any numerical errors which might arise in 
the course of computations v 

A comparison between the structural response of a finite beam and that 
of an infinite one, other physical data being identical, is shown in Fig. 21. 
Only the response near and within the first free-wave propagation zone is 
shown in the figure. Again, the response contains multiple peaks clustered 
roughly within this free-wave propagation zone. However, the number of 
peaks as shown in Fig. 21 exceeds the expected number, that is the number 
of spans of the finite beam. It is believed that some of those peaks have 
resulted from the numerical errors which might have occurred due to the 
ill-conditioned matrices already mentioned. 

(ii) The Two-Dimensional Case 

In this case, the panel is assumed to be made of mylar of which the 
physical properties have been previously mentioned in connection with the 
one-dimensional case. The wires are assumed to be made of steel which has 
the following properties: 


96 




0.736IH 



Fig. 21 Response of a Periodic Beam (Transversely Elastic Supports) 
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10 2 
G = 7.92 3 ; 10^ Wm 

p = 7.68 X 10 ^ kg/m^ 

w 

a =9.65 X 10^ N/m^ 

a 

F = 2.068 X 10^^ N/m^ 

W 


Again, a parametric study using the above values has been carried out 
to obtain an optimal design, in terms of the following physical quantities, 
which meets the requirements set forth in reference [5]: 

(1) the periodic length S, 

(2) the panel thickness h 

(3) the panel width b 

(4) the wire diameter D 

(5) the cavity depth d 

(6) the tension applied to each wire, T. 

Before presenting the results of our parametric study, some introductory 
remarks might be helpful: 

(a) In light of the results of the parametric study carried out for 

“3 

the one-dimensional case, we have chosen I as 2.54 x 10 m. This vajue of £ 
is expected to be near the "best" choice for meeting the wave length 
requirement and, at the same time, keeping the peak region within the 
desired frequency range. 

(b) The wire diameter has been chosen as proportional to the panel 
thickness h, namely 

D = Yh (7.8) 


The proportionality constant y can be chosen arbitrarily but one should keep 
in mind that a large y would make the wires act as rigid supports (a reason- 
able choice might be 2-6). 
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(c) To make sure that the allowable stress of the wire material, a, 

a 

would not be exceeded, the tension T has been chosen as 

(7.9) 

(d) It is to be noted that all computations for structural response 
have been carried out at the midpoint (x = V2, y = b/2) of the first panel. 
Also, the number of modes has been chosen as 5 in the x-direction and 14 in 
the y-di recti on. 

(e) For simpler representation of the numerical results, let us 
rewrite Eq. (5.4) as 


«>^(x, y, (d) = 


H(x, y. 



Sp(aj) 


(7.10) 


/s. 

where H is the "total" frequency response function of the panel and is 
given by 


H(x, y, im ) 


2 


1 

S=-c 


H5(x, 


0) 



(7.11) 


The results of our parametric study are summarized as follows: 

(1) A typical relationship between [Hi and the frequency f is shown 
in Fig. 22 (peaks are not drawn to scale). As can be seen, there are two 
peak regions: the first one occurs at a very low frequency range while the 

second one occurs at a much higher range. When the values (in meters), 
h = 8.98x10“^, I = 2.54x10"^, b = 1.778x10”^ and d ?= 2.54x10”^ were used in the 
computations, the first peak region was found around 10 Hz and the second one 
near 990 Hz. 
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To explain the occurrence of these two peak regions, consider a simpler 
case of a panel of length A and width b which is simply supported along all 
sides and vibrating in vacuum. Assuming a sine series solution in both 
directions, the natural frequency for the nth mode in the x-di recti on and 
the mth mode in the y-direction for such a panel is given by 



where Pp is the panel mass per unit area (= mp/h). Using the values for 
h. A, b and d mentioned above in Eq. (7.12) gives 

fn m = 9-12 {lOOn^ + (m^/49)} 

The first few natural frequencies (i.e., n = 1) are all very close to 912 Hz. 
Therefore, the second peak region in Fig. 22 is mainly due to the panel 
resonance. The presence of the external flow and the cavity appears to 
result in' only minor shift in resonance frequency. 

Secondly, if we consider a simply supported wire of length b, area 
and subjected to a tensile force T, then the natural frequency of the mth 
mode can be computed from 


, -JL / " . ■« 

Pw\ 



(7.13) 


For a wire diameter of .001 inches, the second term under the square root 
sign of Eq. (7.13) is much smaller than the first one. Therefore, 




w 


< 


— /■ ■'^8 = 3 15 m 

2b/ P„ 
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For 1 £m £ 14, the natural frequencies have a range 3-45. Hz. This means 
that the first peak region in Fig, 22 is mainly due to the wire resonance. 

The first and second peak regions of the panel response are plotted in 
much larger scales in Figs. 23 and 24, respectively. 

Since drag reduction is expected to depend primarily on the second peak 
region response, our attention will be focused henceforth on this region and 
any reference to the response will have the restricted meaning of the second 
peak region. 

(2) The panel thickness has a great effect on both the amplitude and 
the location of the peak region. An increase in the thickness h shifts the 
peak region to a higher frequency range and generally reduces the amplitude. 

(3) The effect of the bending rigidity of the wire on the panel 

response represented by the term El — in Eq. (4.36) is negligible. 

w w 

(4) Neglecting the contribution of the torsional moments acting on 
the wire, represented by the term r(m', n, n')in Eq. (4.40), to the virtual 
work slightly reduces the computed response amplitude. 

(5) The change in the tension T has very little effect on the 
structural response of this region (2nd peak), as long as the tensile 
stress in the wire remains unchanged (Thus, the cross-sectional area of 
the wire must be changed correspondingly). 

(6) The change in the panel width (1.778 'x- 3.048) x lo'^m significantly 
alters the locations and the amplitudes of the individual peaks of this 
peak region (see Figs. 24 and 25). 

' (7) No significant change in the structural response was noticed when 

the cavity depth was changed (2.54 i.27)xio~^m (see Figs. 24 and 26). 

Again, it is to be noted that the following values have been used 
in computing the response shown in Figs. 22-26: 
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£ = 2,54 X 10”'^in 
h = 8.98 X 10 ^ m 
j = 4.448 X 10”^ N 
Q = 2 . 54 X 10 ^ m 

I 

The root mean square value of the displacement, _ _ , is computed 
from Eq. (5.6 ). A typical W of 1.74 x io“^ m has been obtained 

I ♦ III • S * 

2 

when using the above data and choosing a cavity depth d = 2.S4 x 10 m and a panel 
width b = i. 778 xiQ^m. This value of W is based on the second peak 

r • III • o • 

frequency range 970-1045 Hz, and it can be reduced by including structural 
damping in the computations. Therefore, damping materials or devices can 
perhaps be used to adjust the response level to within the desired range. 

Finally, the spectral density of the induced pressure P-j is 

plotted in Fig. 27 over a frequency range of 5-2000 Hz, and again in 
Fig. 28 in a larger scale within the second peak region. 
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VIII. CONCLUSIONS 

Based on this investigation, we can conclude that: 

(1) It is possible to design a periodic structure which, when excited 
by turbulent boundary layer pressurie' fluctuations, would have a narrow-band 
response. Within this narrow-band 'region (peak region), the response can 
be adjusted to satisfy predetermined requirements by changing some 
structural or flow parameters. Such requirements can be the amplitude, 
propagation speed, and wave length of the dominant component of a packet 

of waves constituting the narrow-band response. 

(2) It has been found from the parametric study that the membrane 
thickness and the periodic length of the structural construction 
significantly affect the location of the peak region, the amplitude, 
and the wave length of the dominant response. Accordingly, by suitable 
choice of these two parameters, the periodic structure would have a 
resonant mode vibrating at a frequency that falls within the energy- 
containing spectral range of the exciting turbulence pressure field, 
leading to a favorable interaction. 

(3) Structural damping of the membrane or other damping devices 
such as viscoelastic foundations can be used to adjust the amplitude of 
the vibration to within a desired range. 

(4) As far as the structural response is concerned, the effect of 
the additional pressure induced by the structural motion above the 
structure is relatively small. However, this conclusion is based on an 
analysis where viscosity in the fluid has been neglected. A more accurate 
representation of the induced pressure can be obtained if the viscous 
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effects can be Included In the governing wave equation. 'At present, 
there are no analytical solutions which take viscous effects into account. 

(5) The wave propagation approach, is most simple and convenient in 
the analysis of Infinitely long periodic structures. However, this 
advantage diminishes greatly for^ finite periodic structures and the 
analysis Is very much similar to that of the transfer matrix method. 



APPENDIX I 


SOME BRIEF COMMENTS ON THE COMPUTER PROGRAMS USED 

In this Appendix, no attempt will be made to reproduce the complete 

i . 

computer programs used or to discuss them in detail. However, it is felt 
useful to comment very briefly on some of the procedures and subroutines 
which have been helpful in obtaining the numerical results. 

(1) In solving the system of simultaneous equations, Eq. (2.33), 

Eq. (4.45) or (4.58), we have used the I.M.S.L. subroutine LEQTIC. 

(2) The eigensystem subroutine package EISPAC was used to obtain 
the eigenvalues of the matrix C" in Eq. (3.54).* 

(3) The solution of the ill-conditioned system of simultaneous 
equations, Eq. (3.66), was made possible by use of the U.O.I. library 
subroutine GAUSZ. This subroutine employes Gauss elimination technique 
(with complete pivoting). 

(4) The inverse of the Vandermonde matrix Q in Eq. (3.53) was obtained 
using the eigensystem subroutine package MINFIT. This subroutine uses a 
singular value decomposition method to solve a system of equations AX = B. 
The inverse of matrix A is obtained as a special case in which B is replaced 
by the identity matrix. 

(5) Integration of Eq. (5.6) was carried out using Bode's 5-point 
integration method [26], namely 

f(x) clx = II (7fQ + 32f^ + 12 f 2 + 32f3 + 7f^) 


* It is possible to solve the generalized eigenvalue problem of 
Eq. (3.50) directly by_the use of the LZ-algorithm [31] (i.e., without 
inverting the matrix Q). 



where 



fj»f(Xo + jh) ; 0.1,2, 
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